Effect of two gaps on the flux lattice internal field distribution: evidence of two length 

scales from fiSH in Mgi_^.Al^B2. 
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We have measured the transverse field muon spin precession in the flux lattice (FL) state of the 
two gap superconductor MgB2 and of the electron doped compounds Mgi_a;Ala;B2 in magnetic fields 
up to 2.8T. We show the effect of the two gaps on the internal field distribution in the FL, from 
which we determine two coherence length parameters and the doping dependence of the London 
penetration depth. This is an independent determination of the complex vortex structure already 
suggested by the STM observation of large vortices in a MgB2 single crystal. Our data agrees 
quantitatively with STM and we thus validate a new phenomenological model for the internal fields. 



The electronic and superconducting properties of 
MgB2 are determined by the presence of two types of car- 
riers (i.e. bands crossed by the Fermi surface), electrons 
in the tt bands characteristic of all graphite-like mate- 
rials, and holes in the nearly two dimensional a bands. 
There is a large consensus on the origin of the large su- 
perconducting Tc = 39 K being due to the high energy 
phonon mode and to the large density of states connected 
with the narrow a bands(2i|- However also tt electrons 
contribute essentially to superconductivity and, in zero 
magnetic field, the transition temperatures of the two 
fluids appear to coincide Q- 

Magnetic properties (e.g. field penetration) of a two 
gap superconductor are not fully established yet. In par- 
ticular the investigation of the field distribution in the 
flux lattice (FL), a specialty of fiSR, is still at an early 
stage in the presence of two gaps, whereas it is theoret- 
ically and experimentally very well developed for single 
gap materialsHli,!!, 011,0. MgB2 is the ideal simple 
material for such an investigation. 

Efforts have been devoted to the intermediate task of 
calculating, Iflj the London penetration depth A within 
the Eliashberg scheme. Another notable handle to- 
ward internal field distribution was offered by STM sin- 
gle crystal measurements with H \\ c. This imag- 
ing technique directly determines a large core dimen- 
sion = 50 nm, at applied fields surprisingly exceed- 
ing the corresponding fJ.oHc2 — 0.13 T. Calculations by 
Usadel equations[l3| justify this observation introducing 
two length scales, fo- and ^tt, within a moderately dirty 
picture of MgB2. Additional evidence for a direct in- 
fluence of two distinct densities of supercarriers on the 
flux lattice and its field dependence was gathered from 



a small angle neutron scattering (SANS) experiment [l^ 
which detects a crossover field iiqH — 0.5 T connected 
with field induced depletion of tt pairs. 

The well characterized alloy Mgi_2;Al2;B2 provides an 
additional experimental handle, since the doping effect 
of the AP+-Mg^+ substitution as well as its infiuence 
on the phonon spectrum, both accurately described by 
band calculations [ij, lla, ll3 lead to the vanishing of Tc 
for X « 0.45. The key feature is that increasing concen- 
trations of Al produce a progressive filling of the a band, 
hence the stoichiometry offers independent experimental 
control over the two superconducting fluids. 

Early ^SR results^] on the pure sample revealed a 
large variation of the second moment cr^ of the muon 
line with external field, up to 0.6 T, attributing the de- 
pendence to low field pinning effects. In extreme type II 
materials cr^ would be proportional to the inverse squared 
London penetration depth. Subsequent studies 
detected a reduction of in Al and C substituted pow- 
ders, consistent with filling of the a band. Recent mea- 
surements |20| extended the field range to 5 T and at- 
tributed (J^{H) to a field dependence of the coherence 
length. We present here extended data on Mgi_a;Al:cB2, 
< X < 0.35 and offer a novel interpretation of the field 
dependence. 

Powder samples of Mg i_^A l^B-;. were prepared accord- 
ing to standard methods l2ll | . The experiments were per- 
formed on loose powders wrapped in thin Al foil as well 
as on a high density (90% of the bulk value) disk pellet 
of pure MgB2. The experiments were performed on the 
GPS and LTF spectrometers of the S/j,S, Paul Scherrer 
Institut, Villigen, Switzerland. 

The samples were cooled in a magnetic field ^qH es- 
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FIG. 1: TF /iSR X ^ data in 0.6 T; a) rotating frame 
asymmetry {u — 65 MHz, symbols), best fit (solid line, Eq.0 
and fit residues (solid line, displaced for clarity); b) FFT of the 
muon asymmetry (symbols) and of the two components of its 
best fit (dash-dot lines, their sum falls under the symbols); 
also shown, but displaced vertically for clarity, are the FL 
lineshape p{-yB) (dashed line, A = 52 nm, ^ = 23 nm) and 
the Gaussian-convoluted FL lineshape (solid line, same A, ^, 
Gaussian half width Aw — 1 MIfz). 



tablishing a flux lattice (FL) below Tc- The muon spin 
precession signal was collected in two detectors, 1 and 2, 
at the opposite ends of an axis perpendicular to H, in 
a transverse field (TF) /iSR setup. The asymmetry in 
the muon decay is obtained from the count rates Ni in 
the two detectors as Ait) = (7Vi(t) - aN2{t))/{Ni{t) + 
aN2{t)); the effective relative count efficiency a is cal- 
ibrated in the normal metallic phase. Figure ^ shows 
a typical precession asymmetry at T = 1.8 K, trans- 
formed, for visualization purposes, to a rotating frame 
at a frequency u = — 5uj, close to the Larmor fre- 
quency ojl = "2,71^ B, where 7 — 135.5 MHz/T is the muon 
magnetogyric ratio and B the modulus of the local field. 
The damping is due to the FL field distribution, since 
muons stopping close to the vortex cores precess faster 
than those midway between two vortices. The standard 
practice 0,01 for polycrystal data is to fit asymmetries to 
Ae~^ cos{Lot + (j)) and to assume the second moment 

B"^ — B^ of the FL field distribution proportional 

(27r7)2AS2. 



AB2 



to the Gaussian second moment, ct^ 

The one-gap calculation^ of the FL field distribu- 
tion p{B) is given by the dashed line in Fig. ^ (where 
V = jB), with its characteristic peak at the saddle point 
between two adjacent fluxons. Also shown is the convo- 
lution of p{B) with a Gaussian, to empirically account 
for FL disorder introduced by vortex pinning and for 
the distribution of demagnetizing fields in the polycrys- 
tal grains. The curves are obtained with reasonable, al- 
though non optimized values of the coherence length, ^, 
and of the London penetration depth, A. 

The fast Fourier transform (FFT) of the precession 
asymmetry (Fig^, filled symbols) reflects^ directly 
p{B) and displays a very similar asymmetric shape. A 



two Gaussian model 

Af,t{t) = Aie"*''"2i/2cos(wii-H0)-t-A2e~*''"'2/2cos(w2i+0); 

(1) 

improves significantly a large fraction of our fits, in par- 
ticular high statistics ones. The FFT of the two terms are 
shown by dashed-dotted lines in Fig. The solid line 
in Fig. ^ actually shows the best fit to Eq. displaying 
negligible fit residuals, A — Afu. The second moment of 
the composite lineshape may be calculated as 

CT^ = aial^+a2(7l2 + i'^i^^i)^i + i'^2-al)ujl-2aia2UJiUJ2 

(2) 

where — Ai/{Ai + A2) are weights and uji the first 
moments of the two components. 

The experimental values of thus obtained may be 
compared with the FL prediction Q for AB^, which is 
derived from the expression of the local field intensity at 
point r 



B{r) ^ f,,HY,Y 



AB^ 



+ g2A2/(l-/i) 

^-q'end-h) 



l+92A2/(l-/l) 



(3) 
(4) 



where ^ and A are divided by Vl — ^, with h = H/Hc2, 
to account 9] for the field dependence of the Ginzburg- 
Landau order parameter. The distribution p(i3) in Fig.^ 
is directly derived from Eq. 13 whose dependence on H 
is essentially due to the reciprocal FL vector q, since the 
FL lattice cell area is equal to <j)i^/ ^qH, where 4>q is the 
flux quantum, hence q oc \/H. 

In the extreme type II case (A ^ ^), since in prac- 
tice H <C Hc2, the second moment 0,0 does not depend 
on cither H or In these conditions the spatial vari- 
ation of B(r)/fioH is governed only by A, and one has 
cr^(/xs~^) = 0.10734A~2(^in~2)^ However, when ^ is just 
moderately smaller than A, the finite size of the flux core 
produces an excluded volume effect '2^ and the result- 
ing decrease of cr^(i?) vs. H is controlled by The 
cutoff function e^'' ^ is a. standard Q approximation 
to the solution of Ginzburg Landau equations 0, which 
involves a modifled Bessel function cutoff. An alternative 
g-V29C approximation was also suggested [23j. 

Figure |2] shows the values of cr^ obtained by Eq. [21 
from T = 1.8 K measurements as a function of fioH, 
for pure MgB2. The high field single Gaussian widths of 
Ref. HO are also included. Their field dependence agrees 
well with ours, despite the different analysis. Inciden- 
tally the difference between values derived from Eq. [21 
and single Gaussian widths is negligible also for our lower 
statistics low field data on pure MgB2. 

The dashed line shows the best fit to Eq.^(xr = 17.5 
per degree of freedom, A = 71 nm, ^ = 9 nm), which is 
poor and fails to reproduce the clear experimental change 
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— Eq.(4) - one gap 




Eq.(6) - two gaps 




— Hao et al. - one gap 
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FIG. 2: nSR second moments in MgB2 from Eq.|5|(«), includ- 
ing the noH > 0.6 T data from R.ef . l2nl (T, single Gaussian 
fit); the dashed, dahed-dotted lines are one-gap best fits (see 
text); the solid line is the two-gap best fit of Eq.El 




FIG. 3: fiSK second moments in Mgi-^, A1^B2 from Eq. |3 
(•, ■, ♦, A, ► respectively for x ^ 0, 0.1, 0.15, 0.2, 0.25, 0.3, 
including data from Ref. l20l (T, x — 0, single Gaussian fit); 
the solid line is the global two-gap fit (see text) of Eq.|S| Inset: 
low field data. 



of slope around 0.5 T (very close to the crossover field 
detected by SANS 13]). The dashed-dotted hne shows 
the effect of a different cut-off function Q , nearly equiv- 
alent to its approximationj23j . which improves the fit 
{Xr = X), but is still poor both at low and at high fields. 
We suggest that the inadequacy of Eq. 0] is due to the 
two distinct densities of pairs evidenced by SANS. 



A detailed calculation ofp{B) like those of Refs. 
has not been worked out to our knowledge for a two-gap 
superconductor. However, the extension of Ginzburg- 
Landau equations to two order parameters ^ indi- 
cates that p{B) ought to be described in terms of a pen- 
etration depth A and two further length scales, a smaller 

and a larger . Each of these parameters corresponds 
to the square root of a diffusion coefficient (the co- 
herence length, for a clean, single gap system). One may 
therefore expect a change of slope in {H) like that ob- 
served, when H exceeds = (/)o/27r/zo^^. 

We therefore propose a simple extension of Eq. O by 
notingPl that, for intermediate fields, B{r) cx A^^ is 
proportional to the superconducting carrier density n. 
Hence it is reasonable to assume 



Bir) 



(5) 



where the weight Wa- — ncy/{na- -|- n^), is proportional 
to the densities of superconducting carriers rio-, /io-.tt = 
H/Hcr^TT, and i?cr,7r = <^o/27r/xo^^ ^. This expression leads 
to a new second moment expression 



AB2 



+(1 - 



(6) 



We notice that Eq 
of previous works 



.^Idoes not follow the same assumption 
(l^ El , where the square root second 
moment oc A"^ was taken as a linear combination of 
two terms. A linear A"^ = ujp ,^w'^ + oj'^ ^wL is indeed 
that provided by two-gap Eliashberg theoryjlflj in terms 
of the superfluid plasma frequencies ujp,ay^Wa- However 
its direct proportionality with cr^ was just assumed in 
analogy with the one-gap case, and this may well be un- 
justified in the multigap superconductor. 

Figure [5] also shows the best fit to Eq. |B1 (solid line, 
Xr ^ 5.8, A = 49(4) nm, w„ = 0.24(4) nm, = 5.1(2) 
nm, = 23(1) nm) which features a change of slope 
around 0.5 T and follows much better the data to high 
fields, as the improved Xr testifies. No futher appreciable 
improvement is obtained by changing the shape of the 
cut-off. 

Our Mgi_a;Ala;B2 data are shown in Fig.|3|. It is ap- 
parent that increasing values of x lead to a reduction of 
(T^ , i.e. to an increase 0, 0| of the penetration depth 
A(a;). This is not surprising since A"^ is proportional to 
an effective superconducting carrier density n, in analogy 
with the simple London model. 

Since the most relevant high field data are missing for 
some Al concentrations, a fully unconstrained fit of all 
four parameter can be performed only for three samples 
{x = 0,0.1,0.3). In these fits and do not vary 
within errorbars. This is expected, since i) the a band 
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FIG. 4: Mgi_a; Ala;B2: a) penetration depth A (•) and length 
Cir (▼) vs. X from the global best fit of Eq. |S| b) effective 
supercarrier density n oc . Lines are guides to the eye. 



is insensitive to impurity scattering and ii) must be a 
function of (tt carriers on their own do not supercon- 
duct). Hence in first approximation and w„ may be 
safely assumed independent of x. 

The best fit shown by solid lines in Fig. corresponds 
to the following strategy: we performed a global best 
fit to Eq. |B1 of the x — 0,0.1,0.3 data, for which high 
field data are available. Thus we optimized A(a;), ^7r(a^) 
and common average values of the other two param- 
eters {wa = 0.29(4),^^ = 5.0(4) nm). Each of the 
X = 0.15, 0.2, 0.25 data set was then fitted independently 
with the same fixed values of Wa, varying only A and 

Our global fit procedure (x^ — 1-8) determines the 
values of A and ^-n shown in Fig.E^L, where A increases and 
f 7r decrease smoothly with expected. The reduction 
of the effective length with doping indicates that the 
presence of Al reduces the tt carrier mean free path. 

Not surprisingly the x = Q value of A is smaller than 
previous determinations [l^ \2{\ . since this parameters de- 
scribes the field profile around a vortex and we are now 
fitting a more complex fluxon structure. Fig. plots 
the effective n vs. x which displays a dramatic drop, ex- 
trapolating to zero already at x = 0.35 (a sample where 
we failed to detect any FL). Note that our data repre- 
sent a powder average, which masks, e.g. a band effec- 
tive mass anisotropy. However, as long as the anisotropy 
{mc/mabY^'^ does not reduce drastically with doping, the 
/iSR lineshape should be mostly sensitive 25] to the com- 
ponent for H. \\ c. 

Our interpretation can be reconciled with that of 
R.ef . [2fil| . whose MgB2 data are included in our fit. The 
authors invoke a field dependence of the (single) coher- 
ence length, which may well actually correspond to a 
field dependent weight of two length scales. However our 
analysis supports explicitly the existence of two scales in 
the system, and it allows a quantitative comparison with 
Usadel calculations, as well as with the STM datajT^. 
As a matter of fact the latter imply a complex vortex 
structure, since they determines a large = 49.6(9) nm 



value, although their measured gap A^. yields = 13 nm 
and the observed upper critical field is consistent with a 

The Usadel calculations [l2| do justify this observa- 
tion. They predict an apparent vortex length scale 
^„ — 2.7^7r + 0.3^cr, which corresponds to 52(4) nm if we 
identify their length scales with our global fit parameters 
(the large uncertainty accounts for the variation of these 
parameters with cut-off and fitting strategy) . The agree- 
ment is very good and it provides quantitative support 
to our model. 

In conclusion we have shown that two length scales 
govern the field distribution inside a two gap supercon- 
ductor, we have measured the average supercarrier den- 
sity as a function of the a band filling and found remark- 
able agreement with Ginzburg-Landau calculations and 
STM determinations of the vortex dimension. This work 
was financially supported by the Federal Office for Ed- 
ucation and Science, Bern, Switzerland. RDR and AP 
acknowledge the support of INFM grant PRA-UMBRA. 
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